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Abstract 

A theory for the structure of isothermal, self-gravitating gas spheres in pres- 
sure equilibrium is developed for softened gravity, assuming an ideal gas equa- 
tion of state. The one-parameter spline softening proposed by Hernquist & 
Katz (1989) is used. We show that the addition of this extra scale-parameter 
implies that the set of equilibrium solutions constitute a one-parameter fam- 
ily, rather than the one and only one isothermal sphere solution for Newtonian 
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gravity. 

We develop a number of approximate, analytical or semi-analytical solu- 
tions which apply in various regions of parameter space. 

For softened gravity, the structure of isothermal spheres is, in general, 
very different from the Newtonian isothermal sphere. For example, as a 
corollary, we demonstrate the perhaps somewhat surprising result that even 
in the complete absence of rotational support it is possible, for any finite 
choice of softening length e and temperature T, to deposit an arbitrarily 
large mass of gas in pressure equilibrium and with a non-singular density 
distribution inside of r for any r > (neglecting effects of changes in the 
equation of state as well as general relativistic effects). 

Furthermore, it is sometimes claimed that the size of the small scale, 
self-gravitating gas structures formed in dissipative Tree-SPH simulations is 
simply of the order the gravitational softening length. We demonstrate, that 
this, in general, is not correct. 

The main purpose of the paper is to compare the theoretical predictions 
of our models with the properties of the small, massive, quasi-isothermal gas 
clumps (r ~ lkpc, M ~ 1O 1O M and T ~ 10 4 fT) which form in numerical 
Tree-SPH simulations of 'passive' galaxy formation of Milky Way sized galax- 
ies (i.e. simulations not involving energy and momentum feedback to the gas 
from supernova explosions, stellar winds, UV radiation from OB stars etc.). 
We find reasonable agreement, despite the neglect of effects of rotational 
support in the models presented in this paper. 

We comment on whether the hydrodynamical resolution is sufficient in 
our numerical simulations of galaxy formation involving highly super-sonic, 
radiative shocks and we give a necessary condition, in the form of a simple 
test, that the hydrodynamical resolution in any such simulations is sufficient. 

Finally we conclude that one should be cautious, when comparing results 
of numerical simulations, involving gravitational softening and hydrodynam- 
ical smoothing, with reality. 
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1 Introduction 



Due to the incredible increase in computing power provided by the computer 
industry over the last decade or so it has become possible to attempt to 
model, by 3-D numerical simulations, the formation and evolution of galax- 
ies using various combinations of gravitational and hydrodynamical codes 
(e.g. Evrard 1988, Hernquist & Katz 1989, HK89 in the following). Since 
the physical problem involves very large density contrasts fully Lagrangian 
codes, like Tree-SPH (HK89) which is completely particle based, are particu- 
larly well suited for this problem. The Tree-SPH code calculates gravitational 
interactions using a hierarchical tree method (Barnes & Hut 1986) and the 
hydrodynamical interactions using the SPH (smoothed particle hydrodynam- 
ics) method (Lucy 1977; Gingold & Monaghan 1977). 

Primarily to suppress effects of two-body gravitational interactions in 
such simulations the gravitational field of a gas or dark matter particles 
is softened, typically by using a Plummer softening kernel (Evrard 1988) 
or a spline softening kernel (HK89) - see section 2 for more details about 
gravitational softening. 

In the 'passive' variant of such simulations (i.e. simulations not involving 
energy and momentum feedback to the interstellar medium due to supernova 
explosions, stellar winds, UV radiation from OB stars etc.) it is generally 
found that several small, massive gas clumps (r ~ lkpc, m ~ 1O 1O M for sim- 
ulations of Milky Way sized galaxies) are formed quite early in the simulations 
and survive for several Gyrs. The high density gas is nearly isothermal with 
a temperature T ~ lO 4 -^. This is due to the form of the radiative cooling 
function used in the simulations. The primordial gas cools quite effectively 
at temperatures T ~ 10 4 — 10 e K, mainly through collisional excitation of H 
and He + by free electrons, but the cooling function is effectively truncated 
below T ~ lO 4 -^ where the density of free electrons rapidly tends to zero 
as the gas recombines. For low gas densities the radiative cooling can be 
suppressed, at the relevant temperatures T ~ 10 4 — 10 6 i^, by up to two 
orders of magnitude by inclusion of effects of a background, ionizing, UVX 
radiation field, but at the high densities characteristic of the small, massive 
gas clumps the cooling function is essentially unchanged (Efstathiou 1992, 
Vedel, Hellsten & Sommer-Larsen 1994, VHSL94 in the following). 

At first glance it seems somewhat surprising that such small, massive, 
isothermal gas clumps can be in quasi-equilibrium at a temperature of T ~ 



3 



10 4 K, since for Newtonian gravity the mass inside of r of a stationary, isother- 
mal, self-gravitating gas sphere is 



m{< r) < r = 3.2 10 7 ^- r M & , (1) 

Hm p G /ii.2 

where ks is Boltzmann's constant, T the temperature, \x is the mean molec- 
ular weight, m p is the proton mass, G is the gravitational constant, T 4 is the 
temperature in units of 10 4 i^, /X1.2 the mean molecular weight in units of 1.2 
and r is in units of kpc in the last expression. Note, though, that, as shown 
by Ebert (1955) and Bonner (1956), the Newtonian isothermal sphere is only 
stable for 



r < 1.72 



hi ; - 2.03 r KING = 1.21 yf'/'i.f »„"•"• kpc , (2) 



\ fim p Gp 



where p is the central density, r KING is given by equation (13) and n is the 
central number density of hydrogen in units of cm -3 . 

The dense, massive gas clumps are generally supported by rotation to 
some extent, but it is obviously of relevance to investigate whether such 
dense, massive gas clumps could be in equilibrium even in the absence of any 
rotational support at all, as a consequence of the gravitational softening. In 
the following it will be demonstrated that indeed they can. 

In section 2 gravitational softening is briefly discussed. In section 3 pres- 
sure equilibrium solutions for isothermal, self-gravitating gas spheres will be 
derived for softened gravity. In section 4 the theoretical predictions will 
be discussed in relation to the results of numerical Tree-SPH simulations of 
galaxy formation. Section 5 constitutes the conclusion and finally in Ap- 
pendix A the softened gravitational potential and gravitational field of an 
infinitely thin, spherical shell is determined. 



2 Gravitational softening 

One popular way of incorporating gravitational softening in numerical, parti- 
cle based simulations is to replace the Newtonian potential, tp(r) = —Gm/r, 
of a point mass of mass m with the spherically symmetric Plummer poten- 
tial fp{r) = —Gmj\Jr 2 + b 2 , where b is the Plummer scale-parameter. From 
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Poissons equation it follows that the point mass is smeared out into a density 
distribution of total mass m 



which is called the Plummer sphere - see e.g. Binney & Tremaine (1987). 
Note that, for b > 0, pp(r) > everywhere. This is to some extent a 
disadvantage, at least if the gravitational interactions are calculated with a 
Tree code - see the discussion in HK89. 

HK89 proposed, inspired by the work of Gingold & Monaghan (1977), to 
soften the gravitational field of a point mass by using the same spline kernel 
as was used in their implementation of SPH. This spherically symmetric 
spline (softening) kernel, originally introduced by Monaghan & Lattanzio 
(1985), has the advantage that it has compact support at r = 2e, where the 
scale-parameter e is denoted the softening length. The point mass of mass 
m is hence smeared out into a density distribution of total mass m 

1 - 1.5m 2 + 0.75m 3 < u < 1 

Ps P line(r) = — 3 \ 0.25 {2-uf 1<U<2 



u > 2 , 



(4) 

where u = r/e. Note that p sp u n e{ r ) — for u > 2, so outside of r = 2e 
the gravitational field and potential are exactly the same as for Newtonian 
gravity. 

In our Tree-SPH simulations we use the HK89 spline softening for the 
reasons given above. We shall consequently adopt this type of gravitational 
softening in the following theoretical considerations also. 



3 Equilibrium solutions for isothermal, self- 
gravitating gas spheres for softened gravity 

We shall restrict the following analysis to spherically symmetric systems only: 
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The force equation for a stationary, isothermal, self-gravitating sphere 
gives 

du VP _ -> 

_ = -_-V„=0 (5) 

where u is the gas velocity, P the gas pressure, p the gas density and the 
gravitational potential. 

Throughout this paper we assume that the hydrodynamics are described 
by the equations of ordinary gas physics. 

In SPH the gas density (and effectively also all other hydrodynamical 
variables like gas pressure, energy, entropy etc.) is estimated on the basis of 
a smoothing kernel, which, for example, can be of the form given in equation 
(4) - see e.g. HK89 and references therein for further details. The character- 
istic smoothing scale is denoted the smoothing length h. When comparing 
the theoretical models, developed in the following, with the small, massive, 
quasi-isothermal gas clumps, formed in our Tree-SPH simulations, we as- 
sume, in this paper, that the theoretical predictions are directly applicable. 
This assumption is valid here, since h is always much smaller than the char- 
acteristic size of clumps. One should note, though, that this assumption will 
not always be valid. 

Assuming an ideal gas equation of state 

P = Nk B T (6) 

where N is the gas number density, equations (5) and (6) yield for Newtonian 
gravity 

^d\np _d<p _ GM(r) 
dr dr r 2 

where 



f = ^ = r'ci , (8) 

pm p 

is a constant, since T in the following is assumed to be constant, and where 
7 is the adiabatic index, c s the sound speed and 

M(r) = 4tt f p{r')r' 2 dr' . (9) 

JO 

In equation (7) we use Newton's first and second theorems, that the 
gravitational acceleration g{r) depends only on the amount of mass inside of 
r . 
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In the case of softened gravity the situation is very different: for a given 
softening length e the softened gravitational acceleration g t (r) depends on 
the density distribution p{r') for r' G [0, r + 2e], so in particular also on the 
density distribution out to two softening lengths outside of r. 

As will be shown in the following, the change of the structure of isothermal 
spheres for softened gravity relative to the Newtonian case is, in general, quite 
dramatic. 

For Newtonian gravity all solutions to equation (7) can be rescaled such 
that equation (7) effectively has one and only one solution in terms of the 
normalized density and radius 

Pking(?) = — , (10) 
Po 

where 

Po = P(r = 0) (11) 

and 

f = (12) 

f KING 

and where 



r KING 



91 0.595 T 4 - 5 p~° 2 5 no ' 5 kpc (13) 



usually is denoted the King radius (e.g. Binney & Tremaine 1987). 

In the case of softened gravity the gravitational field depends on one 
more parameter than for Newtonian gravity: e, the softening length. One 
would consequently expect that all solutions to the force equation (7) can be 
rescaled in terms of a one-parameter family of solutions and, as we show in 
the following, this indeed is the case. 

3.1 Scaling properties of the solutions 

Let (e, T, p Q ) be given and let p(r) be a solution to 

T ^ = ~9e(r) , (14) 
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where 

rr+2e 

g e ( r ) = 4 :7T g f (r,r')p(r')r' 2 dr' (15) 

Jo 

and where g e (r,r') is the softened gravitational acceleration at r due to a 
unit mass, infinitely thin, spherical shell located at r' . In Appendix A an 
expression for g e (r,r') is derived. 
It is straightforward to show that 

g e i(ar,ar') = a~ 2 g e (r,r') , (16) 

where 

a = - . (17) 

Using this it is easy to show that, for parameters (e' = ae,T,p' Q = a~ 2 p ), 

p\r) = a- 2 p{pr x r) (18) 

is a solution to equation (14). In other words if we know all solutions to equa- 
tion (14) for a given softening length e, then all solutions, for any softening 
length e' = cue, can be obtained by a mere scaling. Now, let again (e, T, p Q ) 
be given and let p(r) be a solution to equation (14). It is straightforward to 
show that, for parameters (e, T' = (3T,p' = (3po), 

p'{r) = (5p{r) (19) 

is a solution to equation (14). Again, if we know all solutions to equation 
(14) for a given value of T, then all solutions for any T" = (5T can be obtained 
by a mere scaling. 

Consequently, in general, the solutions to equation (14) can be expressed 
as a one-parameter family in terms of a normalized density 

p(rVo) = - , (20) 
Po 

where 

r' = cTV (21) 

and 

p' = a 2 r l Po ■ (22) 
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3.2 The linear approximation and self-similarity 

Fig. 1 shows the softened gravitational field, g e {^)i for shells located at 
Usheii = 0.0,0.10,0.30,0.50 and 1.00, where u sheU = r shM /e. For u shM < 0.5, 
the gravitational acceleration at r < e/2 due to such a shell is approximately 
linear in r: 

g e (r) ~ -krM + 0(r 3 ), for r < e/2 and r she// < e/2 , (23) 
where it follows from the Appendix of HK89 that the constant k is given by 

Equation (14) can then be solved analytically assuming that the scale of the 
system is less than or comparable to e/2: 

p(r) = p exp(-ar 2 ) , for r < e/2 , (25) 

where 

kM 

a = . 26 
2T 

If, at r ~ e/2, p{r) « po, then from the consistency relation for the total 
mass 

M ~ Mqo = 47rp / exp(-ar' 2 )r' 2 dr' , (27) 
Jo 

it follows that 

M ~ 27rp r(^)a- 3/2 . (28) 
M can then be obtained as 

M * M TAC = (^gV 5 Po /5 = 3.88 10 8 ?T p~ x r < 4 4' 8 M , (29) 

where n is the number density of hydrogen in units of cm -3 , n = n{r = 0) 
and e 3 is the softening length in units of 3 kpc. 

From equations (24), (26) and (29) it follows that 

a = (-) 2 / 5 tt 3 / 5 G 2/5 T- 2 / 5 e" 6 / 5 p 2/5 . (30) 
3 
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We now introduce the characteristic length scale ttac'- 

r TAC = 3 {2a)- 1 ' 2 = 2.74 T 4 ' 2 ^ n - 2 kpc . (31) 
Equation (25) can then be rewritten as 

Ptac{t) = Poex P (-( 2(W3)2 )) , (32) 

or as 

Ptac{t) = - = exp(-4.5f 2 ) , (33) 
Po 

where f = t/ttac- Hence the density distribution is gaussian in r with 
dispersion ttac/3- So for ttac << e & U solutions can be rescaled in terms 
of one and only one solution, as is the case for the Newtonian isothermal 
sphere, though the nature of the two solutions is very different. 

It is sometimes claimed that the size of the small scale, self-gravitating 
gas structures formed in dissipative Tree-SPH simulations is simply of the 
order the gravitational softening length. From equation (31) it follows that, 
in general, this is not correct. 

3.3 The point mass approximation 

From Fig. 1 it can be seen that, for a given softening length e and a shell 
of mass m located at r s h e ii ~ e/2, a point mass of mass m at r' = re- 
sults in a softened gravitational field, which is a quite good approximation to 
softened gravitational field of the shell and for r > e/2 is a significantly bet- 
ter approximation, than the linear approximation discussed in the previous 
section. 

Rewriting equation (14) as 

f d]np = _ M d^ 
dr dr 

where ip pm is the specific potential corresponding to the softened gravitational 
field of a point mass of unit mass, it follows by integration of equation (34) 
that Ppmij') can be expressed as 

M 

P P m{r) = p exp(-((p pm (r) - (p pm {r = 0))-^-) . (35) 
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In the Appendix of HK89 (p pm (r) is expressed as 

W» = ~Gf(r) ; (36) 

the function /(r) is also given as equation (A3) in Appendix A of this paper. 
Equation (34) can hence be rewritten as 

P P m = f = exp((/(r) - f(r = 0)) S£) . (37) 

The consistency equation for the total mass can not be used in this case, 
since it can be shown that 

rr 

M(r) = 47rpo / P pm (r')r' 2 dr' -> oo as r -> oo . (38) 

However, by setting M = M^c in equation (37), one obtains, for systems 
characterized by ttac ~ e, a generally much better approximation to the true 
solution than for the linear approximation - see the following subsection. 



3.4 The general solution 

A numerical, iterative algorithm was developed to obtain exact solutions of 
equation (14). 

The general solutions can conveniently be expressed as 

p(u; u TAC ) = — , (39) 
Po 

where u = r/e and 

u TAC =^f = 0.913 T 4 - 2 ^ ej - 4 n^ 2 . (40) 

In Fig. 2 the solutions are shown, as functions of u, for a large range of 
the utac parameter. As utac the solutions, as functions of u, converge 
towards a vertical line at r = and as Utac —> oo the solutions, as functions 
of u, converge towards a horizontal line at p — 1. 

As u T Ac ~^ 0, then t T ac becomes the characteristic linear scale of the 
isothermal spheres and the solutions converge as 

p(u;u T Ac) -> Ptac{u/utac) , u T ac -> . (41) 
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Conversely, as utac — ¥ oo tking becomes the characteristic linear scale 
of the isothermal spheres and the solutions converge as 

p(u;u TAC ) -> Pking(u/u king ) , u TAC -> oo , (42) 

where uking = tking/z and there is a one-to-one relation between uking 
and mtac : 

ua™g = 0.249 u 2 TAC . (43) 

In Fig. 3 p/pTAC is shown, as function of r/rxAC, for various values of 
utac and, as expected, p/pTAc ^ 1 as utac —* 0. 

In Fig. 4 p/Pking is shown, as function of r/rxiNG, for various values of 
utac and, likewise as expected, p/pking — * 1 as utac —> oo. 

In Fig. 5 p/ppm is shown, as a function of r/r T AC, for various values 
of utac and, as can be seen from the figure, p/p pm — > 1 as utac ~^ 
- as expected - and it is obvious that the point mass approximation, in 
general, is much better than the linear approximation (Fig. 3), as discussed 
in subsection 3.3. 

Fig. 6 shows M/Mtac as a function of u = r/e for various values of 
utac- As can be seen from the figure, M ~ Mtac for utac ~ 0.5 whereas 
for larger values of utac, M becomes increasingly larger than Mtac and 
increases steadily with u. 

This behaviour is expected, since for utac ~ 0.5 the linear approximation 
breaks down and for the point mass approximation, for the general solution 
and for the Newtonian isothermal sphere it can be shown that M{u) — > oo 
as u — > oo. 

Fig. 7 shows, for e = 3 kpc, a typical value used in current Tree-SPH 
simulations of galaxy formation, and (T4///1.2) = 1, the density profiles of 
the general solution for log(no) = -2, -1, .... , 6. Similarly Fig. 8 shows, for 
the same parameters, M(r) in units of 10 8 M Q . 

As can be seen, the masses are in general much larger, at r ~ 1 kpc, 
which is the typical size of the massive, cold gas clumps formed in our 'pas- 
sive' simulations, than would be expected on the basis of isothermal sphere 
solutions for Newtonian gravity - see equation (1). 

As a corollary it follows from equations (29) and (31) that, for given 
physical parameters (e, T), t T ac ~^ and M T ac ~^ 00 as n — > 00. So we 
obtain the perhaps somewhat surprising result that for any e > and any 
positive value of (T4///1.2) - like ~ 1 for typical small, massive, cold clumps 
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formed in Tree-SPH simulations of 'passive' galaxy formation - it is possi- 
ble to deposit an arbitrarily large mass of gas in pressure equilibrium and 
with a non-singular density distribution inside of r , for any r Q > (neglect- 
ing effects of changes in the equation of state as well as general relativistic 
effects). 

4 Comparing with the results of numerical 
Tree-SPH simulations 

To compare the theoretical predictions obtained in the previous section with 
what is actually found in numerical Tree-SPH simulations of 'passive' galaxy 
formation, we performed a series of simulations, starting from vacuum bound- 
ary, top-hat initial conditions, similar to those described in VHSL94. These 
simulations are not very realistic with regards to galaxy formation, but this 
is of no consequence for the present purpose. 

We carried out four simulations. The first three only differed by having 
£sph = 1-5, 3.0 and 6.0 kpc respectively. The fourth was a higher reso- 
lution simulation - it was identical to the second simulation except that 8 
times more gas particles (each with one eighth of the original gas particle 
mass) were used (but note that no extra phases were added to the initial 
conditions). The purpose of this simulation was to check whether the hydro- 
dynamical resolution was sufficient in our simulations. Other details about 
the simulations are given in Table 1. 

In each of the four simulations we selected the three most massive, cold 
gas clumps ((X4//Z1.2) — 1) at time t = 3.2 Gyr. At this time all the cold 
gas clumps were self-gravitating. They typically form at t ^ 1 Gyr in 'mini' 
dark matter halos, but the dark matter is stripped off during the violent 
relaxation phase at first recollapse of the proto-galaxy at t ~ 2.2 Gyr. 

All the cold gas clumps had, at t = 3.2 Gyr, a radial extent of less than 
or order e. The cumulative mass distribution, at t = 3.2 Gyr, for the largest, 
cold gas clump in simulations #1-3, is shown in Fig. 9 and the cumulative 
mass distribution, at t = 3.2 Gyr, for simulations #2 and #4 (low versus 
high hydrodynamical resolution), is shown in Fig. 10 for all three clumps. 

For each cold gas clump we measured the mass M i, s and the central 
density no and determined the mass M teo , which would be expected for an 
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isothermal sphere of this central density and parameters 63 and (T4///1.2) as 
in the simulations (all gas clumps were in the utac ~ 1 regime). The ratio 
M teo /M obs for the individual clumps is plotted in Fig. 11. As can be seen 
from the figure, M teo /M obs is generally well below unity, indicating that ef- 
fects of rotational support, which are not included in our theoretical models, 
is of importance. For all the clumps M teo jM obs is observed to increase as e 
increases. This is to be expected since, for clumps of similar mass and spe- 
cific angular momentum, the ratio q of pressure force to total gravitational 
force scales approximately as e 3 / 2 when q is well below unity and utac ~ 1- 
Furthermore, as the clump masses decrease, M te0 /M f, s is observed to in- 
crease for a given e. This is also to be expected, though on rather more 
complicated grounds: Theoretically one would expect that the average di- 
mensionless spin-parameter < A > depends only weakly on the magnitude of 
the initial density fluctuations later collapsing and forming the 'mini' dark 
matter halos in which the gas clumps subsequently form dissipatively (e.g. 
Barnes & Efstathiou 1987). If one furthermore assumes that the mass of a 
gas clump depends on the circular speed v c of its 'mini' dark matter halo like 
Mocw" a ~ 3, from e.g. the Tully-Fisher relation, then one can show that, 
for a given softening length, q scales approximately as M~ 2 for q well below 
unity and utac ~ 1- 

For the smallest clump in simulations #1-4 the effect of pressure forces 
appears to be comparable to the effect of rotational support. In fact, for 
softened gravity, the relationship between the physical structure of self- 
gravitating, rotationally supported systems and the gravitational softening 
length e is quite similar to what is found in this paper for isothermal spheres, 
as shown in Sommer-Larsen, Vedel & Hellsten (1997). 

Finally a note on hydrodynamical resolution in Tree-SPH simulations: 
From Figs. 10 and 11 it seems - comparing simulations #2 and #4 - that 
the hydrodynamical resolution is sufficient. At later times, when the clumps 
have merged and super-sonic, radiative shocks have occurred, the results of 
simulations #2 and #4 are still very similar - as an example the surface 
density in the resulting gas disks, at t — 5.2 Gyr, is shown in Fig. 12 for 
simulations #2 and #4. This is quite reassuring, because, whereas in nonra- 
diative shocks the density contrast always is less than the Rankine-Hugoniot 
limit of (^-f), where 7 is the adiabatic index, significantly larger density con- 
trasts most probably occur in the radiative shocks and this could potentially 
cause spurious effects, because of effects of limited shock resolution due to 
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the smoothing, inherent in SPH. 

We propose that one should always test a SPH simulation in the manner 
described above for simulations #2 and #4. Such a simple test constitutes 
a necessary (but not, in general, sufficient) condition that the hydrodynami- 
cal resolution is sufficiently high to adequately resolve super-sonic, radiative 
shocks. 

5 Conclusion 

We have developed a theory for the structure of isothermal, self-gravitating 
gas spheres in pressure equilibrium for softened gravity, based on the spline 
softening kernel proposed by Hernquist & Katz (1989). 

Because the gravitational force depends on an extra scale-parameter, the 
softening length e, relative to Newtonian gravity, the solutions constitute a 
one-parameter family, rather than the one and only one isothermal sphere 
solution for Newtonian gravity. 

For softened gravity the structure of isothermal spheres is, in general, 
very different from the Newtonian, isothermal sphere. 

For example we find, as a corollary, the perhaps somewhat surprising 
result, that, for any finite softening length e and temperature T, it is possible 
to deposit an arbitrarily large mass of gas, in pressure equilibrium and with 
a non-singular density distribution, inside of r , for any r > (neglecting 
effects of changes in the equation of state and general relativistic effects). 

Furthermore, it is sometimes claimed that the size of the small scale, 
self-gravitating gas structures formed in dissipative Tree-SPH simulations is 
simply of the order the gravitational softening length. We demonstrate, that 
this, in general, is not correct. 

The main purpose of the paper is to compare the theoretical predictions 
of our models with the properties of the small, massive, quasi-isothermal 
gas clumps formed in numerical Tree-SPH simulations of 'passive' galaxy 
formation of Milky Way sized galaxies (i.e. simulations not involving energy 
and momentum feedback to the gas from supernova explosions, stellar winds, 
UV radiation from OB stars etc.). We find reasonable agreement, despite the 
neglect of effects of rotational support in our theoretical models. 

We would expect that if the gravitational softening was based on a Plum- 
mer kernel (which is also a one scale-parameter type of softening), rather than 
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the spline kernel, the resulting isothermal spheres would be qualitatively sim- 
ilar to those described in this paper. We encourage any Plummer softening 
'fan' to check this quantitatively by going through a similar exercise as de- 
scribed in this paper for the spline softening. 

We have not discussed the stability of the isothermal sphere equilibrium 
solutions obtained in this work - we shall return to this issue in a forthcoming 
paper using an approach similar to that outlined by Ebert (1955) and Bonner 
(1956). 

We comment on whether the hydrodynamical resolution is sufficient in 
our numerical simulations of galaxy formation involving highly super-sonic, 
radiative shocks and we give a necessary condition, in the form of a simple 
test, that the hydrodynamical resolution in any such simulations is sufficient. 

The results obtained in this paper, and in Sommer-Larsen, Vedel & Hell- 
sten (1997), seem to indicate that one should be cautious when comparing 
what is observed in the real Universe with results of numerical Tree-SPH sim- 
ulations or any other numerical, gravitational-hydrodynamical simulations, 
where the calculation of the gravitational interactions, locally or globally, is 
particle based and the gravitational field of the individual particles is soft- 
ened. 

To summarize: Without proper testing/understanding of the effects of 
gravitational softening and hydrodynamical smoothing in numerical gravitational- 
hydrodynamical simulations of various physical problems, the degree of real- 
ity represented by such simulations may be quite difficult to assess. 
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Appendix A: The softened gravitational po- 
tential and gravitational field of an infinitely 
thin, spherical shell 

It is fairly easy to show that the potential, at radial coordinate r, of an 
infinitely thin, spherical shell with radius R and unit mass is given by 

2k fg% <P P m(r'ydr' , r<R 

<Psheii(r,R) = I (Al) 

m Ir T -R V P m(r')r'dr' , r>R 

where <p pm is the softened potential of a point mass of unit mass. In the 
Appendix of HK89 (p pm is expressed as 

V P m(r) = -Gf(r) , (A2) 

where 

-[2u 2 /3 - 3w 4 /10 + w 5 /10]/e + 1.4/e < u < 1 

fir) = | -l/(15r) - [4M 2 /3-M 3 + 3M 4 /10-M 5 /30]/e + 1.6/e \<u<2 

1/r u>2 

(A3) 

The gravitational acceleration, at r, from the unit mass shell at R, is then 
given by 



g e (r,R)= - dip she u /dr 



(A4) 



and 



ihf£-r r <P P m(r')r>dr> 
-^(v P m(R + r) - (f pm {R - r)) 
-^{ippm{R + r) + ip pm (R - r)) 



r < R , 
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lh Ir-R^pm{r')r'dr' 

-^(Vpmir + R) - ippm(r - R)) 

-^(<Ppm(r + R) + fpm(r + R)) r>R. 

It is possible to calculate g e (r,R) analytically, but as the result is quite 
complicated and voluminous, we shall omit it here. 

From equation (A4) it follows that for a spherical system, with mass 
distribution p(r), the gravitational acceleration at r is given by 

g e (r) = 4tt f ^ g e (r, r')p(r')r' 2 dr' . (A5) 
Jo 
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Figure captions 



Figure 1: Gravitational acceleration due to infinitely thin, spherical shells of 
mass m located at u s h e ii = T s h e ii/t = 0.00, 0.10, 0.33, 0.50 and 1.00. Solid line 
corresponds to u s h e ii = 0.00 and as u s h e ii increases the curves move monotoni- 
cally upwards. The linear approximation is shown as the straight dotted line. 

Figure 2: The general solutions for isothermal spheres for softened grav- 
ity. As utac increases the solutions become monotonically more horizontal. 

Figure 3: p/ptac as a function of r/r T AC- The solid, horizontal, straight 
line corresponds to utac — 0.00. As utac increases the curves move mono- 
tonically towards the other solid line, which represents utac = oo. 

Figure 4: p/pking as a function of t/tking- As utac increases the curves 
converge towards the solid, horizontal, straight line, which corresponds to 
utac = oo. The other solid line corresponds to utac = 0.00. 

Figure 5: p/p pm as a function of t/ttac- As utac decreases the curves 
converge towards the solid, horizontal, straight line, which corresponds to 
utac = 0.00 

Figure 6: M/M T ac as a function of r/e. As utac increases the curves start 
out monotonically less steep. 

Figure 7: Density distributions for isothermal spheres, of various central hy- 
drogen number density no, for T = 10 4 K and gravitational softening length 
e = 3.0 kpc. 

Figure 8: Cumulative mass distribution for isothermal spheres, of various 
central hydrogen number density n , for T = 10 4 K and gravitational soft- 
ening length e = 3.0 kpc. The lines correspond to the same values of n as 
in Fig. 7. 

Figure 9: Cumulative mass distribution, at t — 3.2 Gyr, for the most massive, 
cold gas clump in simulations #1-3. 
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Figure 10: Cumulative mass distribution, at t — 3.2 Gyr, for the three most 
massive, cold gas clumps in simulations #2 and #4. 

Figure 11: The ratio of the theoretical to observed mass, at t — 3.2 Gyr, 
for the three most massive, cold gas clumps in simulations #1-4. For more 
details - see text. 

Figure 12: The gas surface density distributions, at t — 5.2 Gyr, for the 
disks formed in simulations #2 and #4. 
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Table 1: Tree-SPH simulations #1-4. 



Simulation e SPH £dm N gP # N dm 

[kpc] [kpc] 

#1 ~ 1.5 10.0 2048 2048 

#2 3.0 10.0 2048 2048 

#3 6.0 10.0 2048 2048 

#4 3.0 10.0 16384 2048 
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